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Summary 

The  current  Interest  In  rockets  and  space  travel  has 
aroused  a  corresponding  Interest  In  the  determination  of 
majclmum  range,  minimum  time,  and  so  on,  for  various  types  of 
trajectories . 

A  variety  of  questions  of  this  type  have  been  treated 
by  means  of  the  theory  of  dynamic  prograjnndng .  Here  we  wish 
to  show  how  to  use  functional  equations  to  determine  the 
range,  the  maximum  elevation,  and  similar  quantities,  as 
functions  of  Initial  position  aiid  velocities. 
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FUHCTIONAL  EQUATIONS  AND  MAXIMUM  RANGE 
Richard  B«llman 


1.  Introduction 

'Rie  current  interest  in  rocketa  and  space  travel  has 
aroused  a  eorx*e spending  interest  in  the  determination  of 
oaxlmun  range,  mininnun  time,  and  so  on,  for  various  types  of 
trajectories. 

A  variety  of  questions  of  this  type  have  been  treated 
by  aeauis  of  the  theory  of  dynamic  piegramming,  see  j^l,2,4|. 
Here  we  wish  to  show  how  to  use  functional  equations  to 
determine  the  range,  the  maxinum  elevation,  emd  similar 
quantities,  as  functions  of  Initial  position  and  velocities. 

2.  Vertical  notion— I 

Consider  an  object,  subject  only  to  the  foirce  of  gravity 
and  the  resistance  of  the  air,  which  is  propelled  straight 
up.  In  order  to  illustrate  the  technique  we  shall  employ, 
let  VIS  treat  the  problem  of  determining  the  maximum  altitude. 
Let  the  defining  equation  be 

(1)  u"  -  -  g  -  h(u' ), 

with  the  initial  conditions  u(0)  -  0,  u'(o)  =  v.  Here 

V  >  0,  and  h(u'  )  ^  0  for  all  u'  . 

Since  the  maximum  altitude  is  a  function  of  v,  let  us 
introduce  the  function 


D 


(2) 
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f(v}  -  the  maxlfflum  altitude  attained  staining 
with  initial  velocity  v. 


Prom  the  definition  of  the  fvmctlon  It  follows  that 

(3)  f(v)  -  vA  +  f(v  -  [^g  +  h(v)]A)  +  o(a), 

for  A  an  Infinitesimal.  Verbally,  this  states  that  the 
maximum  altitude  is  the  altitude  gained  over  an  Initial 
time  A,  plus  the  maximvun  altitude  attained  starting  with  a 
velocity  V  —  fg  +  h(v)jA,  the  velocity  of  the  object  at 
the  end  of  time  A,  to  within  o(a). 

Expanding  both  sides  and  letting  A  — >  0,  we  see  that 

(4)  f.(v)  =  - 1—  . 

g  +  h(v) 


Since  f(0)  -  0,  this  yields 


(5) 


f(v)  ./ 


V  'l'*''! 


0  g  +  h(vj^) 


In  the  particular  case  where  h(v)  =  0,  we  obtain  the 

2 

stemdard  result  v  /2g. 


3.  Vertical  Motion--II 

Consider  the  more  general  case  where  motion  la  through 
an  Inhomogeneous  medium.  Let  the  defining  equation  be 


'2* 


(1) 


h(u,u'),  u(0) 


u'(0) 
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Absxum  that  h(u,u' )  ^  0  for  all  u  aiid  u',  so  that 
-  0  Implies  no  motion. 

The  maximum  altitude  Is  now  a  function  of  both  and 

C2>  Introduce 

(2)  f(c^,C2)  »  the  maximum  altitude  attained  8tax*tlng 

with  the  Initial  position  c^^  and 
Initial  velocity  c^. 

Then,  as  above, 

(3)  f(o^,C2)  =.  C2A  +  f(c^  + 

which  yields  In  the  limit  the  partial  differential  equation 

®2  ^  °2  h(c^,C2)  “  0- 

By  virtue  of  our  assumptions,  f(cj^,0)  s  0,  for  c^^  ^  0. 

4.  Computational  Aspects 

One  can,  of  course,  use  the  method  of  characteristics, 
or  standard  difference  methods,  to  solve  (3*4).  Let  us 
present  another  method  which  reduces  the  solution  to  the 
tabulation  of  a  sequence  of  functions  of  one  variable. 

In  place  of  (3-4),  let  us  use  the  discrete  approxl- 
BULtlon  of  (3*3 ), 

(1)  f(o^,C2)  =02^  +  f(c^  +  °2^'®2 
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Slno*  iB  Bionoton*  d«or«aslng.  It  oan  4«  us«d  to  play 

the  role  of  time.  Let  ue  write  n6, 

where  b  la  a  positive  quantity,  and  f(cj^,02)  «  fj|(Oj^). 

We  consider  then  only  values  of  Cg  which  are  multiples  of 
b.  To  overcome  the  fact  that  Cg  +  h(cj^,02)A  In  general 
will  not  be  a  multiple  of  b,  we  can  either  replace  It  by 
f(cg  +  h(c^,Cg)A )/b] ,  or  use  Interpolation.  Although  use 
of  an  Interpolation  formula  slows  up  the  computation.  It 
greatly  Improves  the  accuracy.  For  an  application  of  the 
foregoing  techniques  to  a  more  complicated  partial 
differential  equation,  see  [s]. 

5.  Maximum  Altitude 

Consider  now  the  case  where  motion  takes  place  In  a 
plsuie.  Let  the  equations  be 

(1)  X"  -  g(x',y'),  x(0)  =  0,  x'(0)  -  Cj^, 

y"  -  h(x',y'),  y(0)  «  0,  y'(0)  -  Cg. 

Introducing,  as  before,  the  function  f(cj^,Cg)  equal  to  the 
maximum  altitude,  we  see  that 

(2)  f(Oj^,Cg)  -  (cj  + 

+  f(Cj^  +  g(Cj^,Cg)A,0g  +  h{cj^,Cg)A)  +  o(a). 


Hence, 


(3) 


(c^  +  +  g(o^,Cg)  +  h(c^,Cg)  ^  -  0. 
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Once  a^aln,  let  us  assur.ie  that  Cg  =•  0  Implies  no  vertical 
motion.  Then  fic^,0)  .  0  for  ^  0.  It  follows  that  we 
can  again  compute  the  solution  by  means  of  a  sequence  of 
functions  of  one  vairlable. 

6 .  Maxlmun  Range 

To  tackle  the  problem  of  maximum  z*ange  directly  requires 
the  introduction  of  another  state  vartable,  the  Initial 
altitude.  It  can  easo  be  broken  up  Into  two  problems, 
corresponding  to  the  ascent  to  maximum  altitude,  and  the 
descent . 
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